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Thermal vector potential formulation is applied to study thermal dynamics of magnetic structures 
in insulating ferromagnets. By separating variables of the magnetic structure and magnons, the 
equation of motion for the structure including spin-transfer effect due to thermal magnons is derived 
in the case of a domain wall and a vortex. The magnon current is evaluated based on a linear response 
theory with respect to the thermal vector potential representing the temperature gradient. It is 
shown that the velocity of a domain wall when driven by thermal magnon has a strong temperature 
dependence unlike the case of an electrically-driven domain wall in metals. 


I. THERMAL EFFECTS IN SPINTRONICS 

Manipulation of magnetization and spin current with¬ 
out magnetic field is a key issue in spintronics technol¬ 
ogy. In the case of metallic ferromagnets, electric cur¬ 
rent has been proved to be highly useful for switching 
the magnetization and for driving magnetic domain walls 
by use of current-induced spin-transfer torques proposed 
by Berger-ii^ and Slonczewski^. Metallic systems have 
been studied extensively also for spin current manipula¬ 
tion such as for direct and inverse spin Hall effects^i^. 
Besides metallic magnets, insulating ferrimagnet such as 
yttrium iron garnet (YIG) are also expected to be useful 
spintronics materials^ because of its weak spin relaxation 
effects as represented by small Gilbert damping param¬ 
eter, a. Although use of electric current do not apply 
to insulators, recent extensive studies revealed that ma¬ 
nipulation of magnetization of YIG is possible by several 
different methods, such as using a temperature gradient^ 
and sound waves^. Thermal methods are of particular 
importance for realizing novel thermoelectric materials 
based on magnetic materials^ii^, which are expected to 
be useful for devices with low energy consumption. 

Recently, thermally-driven domain wall motion in a 
YIG film was observed opticallyii. The domain wall was 
found to move to a hotter side of the system, and the 
speed was about 2 x 10“^ m/s when a temperature gra¬ 
dient of VT = 20 K/mm is applied. The direction of 
the motion is counter-intuitive, but is consistent with the 
spin-transfer effect of spin waves (magnons) ^^’^^ . In fact, 
magnons excited around the wall are pushed by the tem¬ 
perature gradient to the colder side, and they transfer 
an spin angular momentum opposite to the local magne¬ 
tization, resulting in a domain wall motion towards the 
hotter end. Thermodynamic viewpoint of domain wall 
motion was discussed by Wangi^ by evaluating numer¬ 
ically the entropy and free energy due to the magnons 
around the wall. It was concluded that the wall motion 
to the hotter side is consistent with minimization of free 
energy when the entropy due to magnons are taken into 
account. In a metallic ferromagnet, thermal motion of 
domain walls was observed in 1986^^. The motion was to¬ 
wards the colder end, and the wall-entropy force propor¬ 


tional to VT was suggested as a possible mechanis m^^d^ , 
but the behavior would be explained simply by the spin- 
transfer torque induced by thermally-driven conduction 
electrons. 

Domain wall under temperature gradient was numer¬ 
ically studied by solving the Landau-Lifshitz-Gilbert 
(LLG) equation for spins in Refii^, where the effect of 
temperature was modeled as a random magnetic field 
satisfying the fluctuation-dissipation theorem. The re¬ 
sults indicated that the behavior of the wall are essen¬ 
tially the same as the electric-current-driven casei^^— . 
For instance, a sliding motion with a constant tilting an¬ 
gle occurs at low driving force or at large damping cases, 
while the tilting angle becomes time-dependent result¬ 
ing in a screw-like motion when the driving force exceeds 
a threshold value corresponding to the Walker’s break 
down field. The LLG approach was employed also to 
study the motion of magnetic skyrmions^i. Skyrmions 
also turned out to move to the hotter region in the same 
manner as domain walls and the velocity was found to be 
inversely proportional to the Gilbert damping parameter, 
a. The numerical result was explained by the “Brownian 
motion” of the skyrmion driven by thermal random mag¬ 
netic field and the spin-transfer effect from thermal spin 
waves. It was demonstrated that the spin wave spin- 
transfer effect contributing to the motion towards the 
hotter side is proportional to VT/a, while the velocity 
due to the Brownian motion, which pushes the skyrmions 
to the colder side, proportional to aVT, and thus is neg¬ 
ligibly small when a is small. The study of Refi^ was 
extended to include dissipative torque arising from spin 
relaxation in Refs i^^’^^ . 


II. LUTTINGER’S FORMULATION OF 
THERMAL TRANSPORT 

Most of the theories on thermally-driven magnetic 
systems so far are either phenomenologica l ^ i ^ or 

numericali^i^ and lack microscopic viewpoint. The rea¬ 
son is obvious; the temperature gradient appears not to 
be straightforwardly integrated into a linear response the¬ 
ory or field-theoretical methods since its effect is not de¬ 
scribed by a microscopic interaction Hamiltonian. This 
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difficulty was removed in the case of electron transport 
phenomena by Luttinger already in 19643^. He intro¬ 
duced a hctitious scalar field called a ‘gravitational’ po¬ 
tential, 'b, which couples to the local energy density, f, 
by an interaction Hamiltonian 

Hl = J (fr'SS. ( 1 ) 

It was argued that 4/ satisfies V'b = ^2^ and that a linear 
response theory (Kubo formula) applies to the thermally- 
driven case by considering the correlation functions of 
the energy current density. Many works have been car¬ 
ried out based on the ‘gravitational’ potential formal¬ 
ism on electron transport33“— , Quantum dol33, magnon 
transport3331 and thermal torques^. At the same time, 
it has been noted that naive application of Kubo formula 
may result in wrong thermal coefficient923313i33. 

As alternative approaches for thermal transport, 
Landauer-like approach by connecting locally equilibrium 
systems at different temperatures is possible^3£. It was 
shown in Ref.— that the obtained transport coefficients 
have no unphysical divergence, since only the electrons 
at the Fermi energy contribute to them. Quantum ki¬ 
netic equation approach was employed in Refs— 33 to 
study electron correlation effects. Still, a fully quantum 
mechanical description based on a linear response theory 
has wider applications and is highly useful. 

Prescriptions to extract correct results based on the 
‘gravitational’ potential formalism were given by some 
worka333i33. In the case of thermally-induced electron 
transport, Smrcka and Streda carried out the calcula¬ 
tion of thermal Hall effect in the presence of an applied 
magnetic field, based on the Luttinger’s prescriptior33i21. 
They rewrote the Luttinger’s interaction, 'kf, by use of 
integral by parts as —£r ■ and carried out a pertur¬ 
bative expansion with respect to a physical force propor¬ 
tional to V'k = VT/T. The analysis, however, is based 
on a interaction proportional to a unbounded operator, r, 
and thus the results may not be always convincing. Re¬ 
cent works on thermal magnons are based on the same 
treatment^!. Later, Qin et. al. pointed out that a naive 
application of Luttinger’s approach for thermal Hall ef¬ 
fect leads to a unphysical divergence at T —^ 0, and that 
this divergence arises from an equilibrium rotational elec¬ 
tron current induced when the time-reversal invariance 
is broken by the magnetic field^i. They showed that the 
correct non-equilibrium response is obtained if one sub¬ 
tract the equilibrium current before applying the ‘gravi¬ 
tational’ potential. 

The case of thermally-induced torque in ferromagnetic 
metals was studied in detail by Kohno et. al— . They 
calculated the non-equilibrium torque on the magnetiza¬ 
tion generated by the conduction electron when a tem¬ 
perature gradient is applied and found that unphysical 
divergence arises if Luttinger’s approach is straightfor¬ 
wardly applied. They demonstrated that the divergence 
arises from the equilibrium torque describing the ex¬ 
change interaction between the magnetization and that 


this equilibrium contribution needs to be removed when 
discussing non-equilibrium torque. 


III. VECTOR POTENTIAL FORMULATION 

It was recently pointed out in Ref— that those prob¬ 
lems of Luttinger’s formalism arise from the fact that 
the ‘gravitational’ potential couples to the total energy 
density, modifying the equilibrium properties in addition 
to inducing non-equilibrium response. It was discussed 
that the role of diamagnetic current, which is essential for 
removing unphysical equilibrium (non-dissipative) con¬ 
tribution from transport coefficients3i33^ is not seen 
straightforwardly in the scalar potential formalism. In¬ 
stead of a scalar potential formalism, a vector potential 
formalism to describe thermally-induced transport was 
developed in Refs— 33. In Ref—, Shitade demonstrated 
using the analogy of general relativity that if an invari¬ 
ance under time translation is imposed locally, a vec¬ 
tor potential arises from Luttinger’s scalar potential and 
they are described by a gauge invariant theory. He ap¬ 
plied his model to describe thermal Hall effect of non¬ 
interacting electrons and showed that the results satisfy 
the Wiedemann-Franz law. The origin of the invariance 
under local time translation was not argued. In Ref—, 
the vector potential was introduced by rewriting the Lut¬ 
tinger’s Hamiltonian for the ‘gravitational’ potential us¬ 
ing the conservation law of the energy in the limit of 
static (dc) limit of the temperature gradient. The vector 
potential interaction was also argued there in the context 
of entropy force. Since the thermal vector potential cou¬ 
ples to energy current, it only generate excitations and 
does not alter equilibrium contributions. The vector po¬ 
tential representation thus guarantees a straightforward 
linear response calculations for thermal dynamics on the 
equal footing as the electric field-driven case. It was in¬ 
deed shown that unphysical equilibrium contributions are 
canceled automatically by ‘diamagnetic’ currents associ¬ 
ated with the vector potential. A possibility of vector 
potential description was briefly mentioned in Ref— in 
discussing magnon-drag thermoelectric effects. 

The vector potential form of the interaction Hamilto¬ 
nian describing the thermal effect is 

Hat = ~ J • Arit) (2) 

where Axit) is the thermal vector potential. It satisfies 

dtArir, t) = V^'(r, t) = (3) 

where tk is the Luttinger’s ‘gravitational’potential. The 
interaction ([2]) was used to describe thermally-induced 
longitudinal transport and Hall effect of non-interacting 
electrons3i33. In the present paper, we apply the formal¬ 
ism to study thermally-induced magnetization dynamics. 
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IV. ENERGY CURRENT OF LOCALIZED 
SPINS 

To study thermal transport based on Eq. 0 , we need 
to derive the expression for the energy current density. 
Following Ref4i, we carry out the derivation quantum 
mechanically by use of energy conservation law, 

£ + V-j£=0, (4) 

where £ is the energy density. We consider a case of a 
ferromagnet with an easy-axis and a hard-axis magnetic 
anisotropy energies. The Hamiltonian is H = J ^£, 
where 

f = + (5) 

where J, K{> 0) and A"_l(> 0) represent the energy of 
the exchange interaction, the easy-axis anisotropy and 
the hard-axis anisotropy, respectively, and a is the atomic 
lattice constant. The energy current density in this case 
is obtained as (see Sec. El for derivation) 

3£,^ = -^y^S■S. ( 6 ) 

Thermally-driven dynamics of the spin structure is calcu¬ 
lated below as a linear response to the interaction H^t i 
E q. dH), with the energy current density of Eq. ([S]). 


where 7 = y ^. We neglect the terms third- and higher- 
order in boson operators. 

The unitary transformation modifies derivatives of spin 
as {fi = t,x,y,z) 

d^S = U{d^+iAjj)S, ( 10 ) 

where 


( 11 ) 

is a spin gauge field represented by a 3 x 3 matrix. Ex¬ 
plicitly, the spin gauge field reads 




0 0 

1 \ 

/ ° 

- COS0 

0 

0 0 

0 

— iVncj) cos 9 

0 

sin 

-1 0 

0/ 

V 0 

— sin 6 ( 

0 




(12) 


The energy current density, Eq. ®, is then 
je,i = — nS{dt —iAu + iAu i)S 
= S^{{v,e)0 + sin^ e{v,4>)^) (^i - 

+ s(bi{V,b) + iVd 


-i cos 9[^{bi V b) + (Vi(())(6'^ dt &)] 


V. SEPARATION OF VARIABLES 

To carry out the calculation of thermally-driven mag¬ 
netization dynamics, we separate a collective degrees 
of freedom describing a classical magnetization struc¬ 
ture and fluctuation, magnons or spin waves, around 
the structure. The direction of the localized spins for 
the classical solution is represented by use of polar an¬ 
gles 0{r,t) and The magnon excitation around 

the structure is then represented by use of the Holstein- 
Primakov boson defined with respect to the local quan¬ 
tization axis along The localized spin 

vector, S, is thus represented as 

S = U{r,t)S = U{r,t){Sz + Ss), (7) 

where z is the unit vector along the z axis, C/ is a 3 x 
3 unitary matrix describing a rotation of a vector z to 
the direction {9, (j)), S = Sz + Ss, and 6s represents the 
fluctuation. The unitary matrix is chosen as^ 

( cos 9 cos (j) — sin ^ sin 9 cos (j) \ 
cos 0 sin(/) coscf) sin^sint/) , ( 8 ) 

— sin 9 0 cos 9 J 

and the fluctuation is represented in terms of annihila¬ 
tion and creation operators for magnon (the Holstein- 
Primakov boson), b and b\ as^ 

/ 7(6'I' +b) \ 

Ss = i7(6t - 6) , (9) 

V -fctfo J 


+2 cos^ 9{Vicj))^b% 

+ Sj \-i cos 6i[0(V,(/)) -I- (V*6»)0] (fel' - b) 

+ {9 + ism9^){'Vib'^) + {9 — isin00)(Vi&) 
-t-(Vi0 -I- i sin 9'Vi(j)){b^) + (V^fl — i sin 9'Vi(j)){b) 


(13) 


The first term, 

jg = -^S^{{W,9)9 + sin 2 0 (V,(^)^) (^1 - ^b^b'^ , 

(14) 


is the energy current carried by the magnetization struc¬ 
ture, and 


•(m) _ 

rs,^ = 


JS 


6 t(V, 6 ) + (V, 6^)6 


(15) 


is the contribution of magnons. We see that the tempera¬ 
ture gradient appears to act on the magnetization struc¬ 
ture as well as on the spin waves. We need to be careful, 
however, since Eq. HMD may contain a steady-state con¬ 
tribution not contributing to the thermally-excited re¬ 
sponse (see Sec. Nm . The mixed contributions in Eq. 
(1131) where the energy current is carried by both the spin 
texture and the magnon turn out to be second order in 
the temperature gradient and are thus neglected. The 
last contributions linear in the magnon operators are ne¬ 
glected similarly. It is thus sufficient to consider the two 
energy currents, Eqs. (HI (US]) in the following analysis 
of linear response regime. 
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VI. MAGNON SPIN-TRANSFER EFFECT 


The unitary transformation gives rise to a gauge field, 
Ajj, due to the exchange interaction according to 

{VSf = iysf - iS\Au- V)S + S^{AufS, (16) 


where V^=V^ — acts on field operators on the 
both sides but not on Au, i.e., S^Au- V)S = 
Aj/The exchange interac¬ 
tion term is thus 


J 

2 


(VS)’ = i 



{\7Sf+Hf, 


(17) 


where 


is ‘diamagnetic’ current of magnon. The total magnon 
spin-transfer effect is thus described by 


Hst = US 


(firAs ■ jm, 


(23) 


where jm = jm ^ + jm^ ■ 


VII. DOMAIN WALL 

Let us describe a domain wall based on the Hamil¬ 
tonian (ID). Considering a thin wire, the magnetization 
structure is treated as one-dimensional, i.e., changing 
only in the wire direction, which we choose as z axis. 
The classical solution is determined by the equation of 
motion as 


Hf ^ J ^SHAu- y)S 

= US J d^rAs-ji^\ (18) 

represents the interaction between the spin structure and 
the magnons (subscript st denotes spin-transfer), and 
contributions second order in Ajj are neglected. Here 
As= \ cos 6\7(j) is the adiabatic component of spin gauge 
field^S and 


V&), (19) 

is the magnon current density. As seen in (TTSl) . the 
magnon current couples to the adiabatic spin gauge field 
in the same manner as the case of the electric current in 
metals in the adiabatic regime^Siil. The magnon current 
therefore induces the spin-transfer torque with the an¬ 
gular momentum transfered per area and per unit time 
is 4J5'ji„ (instead of in the case of electric current 
{P is the spin polarization of the current and e is the 
electron charge)). 

There is another contribution to the interaction be¬ 
tween the magnetization structure and magnons, namely, 
the one arising from Hat (Eq.(I2|)). Using Eq. (IT^ . the 
spin-transfer terms in Hat is 


cos0 = tanh^^—sin0 =- —wtt, (24) 

A cosh^^ 

and Vz4>{t) = 0, where A = is the thickness of the 
domain wall. We then have 

WzS = -jee, (25) 

where eg = (cos 6 cos (/), cos 9 sin </>, — sin 9). 

We first show that the energy current of the wall given 
by Eq. (1141) is an equilibrium contribution, which should 
not be taken into account. Let us first see what happens 
if Eq. (IT4l) is naively applied. Equation (IT4l) for domain 
wall reads (neglecting the corrections of the order of S~^ 
by magnons) 


'^*’"A2a3 cosh^i^' 


(26) 


The coupling between the domain wall and temperature 
gradient, Eq. then becomes 


Ha^ = 


2AJS‘^ VzT 


dt'X{t') = V, 


JS^ VzT 

■^2 Y' 


A, 

(27) 


H% = i^AT,r J d\[{dtS)Au,,S + {X,S)AuAS]. 

( 20 ) 

Keeping only the contributions second order in the 
magnon operators, it reads 

Ht=UsJd^rAs-ji^\ ( 21 ) 

where 

^ -^Arifi K &), ( 22 ) 


where = 2XA/Y is the number of spins in the wall 
{A is the cross sectional area of the wire). One may be 
tempted to conclude that the thermal force on the wall 
is 


Fx = 


5Ha^ ,, v,r 

5X "2hX^ T 


(28) 


As we see, this force diverges at T —0 if VT is fixed, 
which is unphysical. The problem arose from the energy 
current density of Eq. (l26l) . Equation (l26l) indicates that 
the energy current arises proportional to the wall speed 
X. This energy current is, however, a steady-state cur¬ 
rent, which does not contribute to the entropy change. 
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since the energy carried by the domain wall is not con¬ 
verted to a heat unless the wall is annihilated. To con¬ 
sider the standard experimental situation where the wall 
is semi macroscopic and does not annihilate, one needs 
to thus subtract the steady-state contribution from Eq. 
()26p . resulting in vanishing energy current for the wall. 
This treatment is consistent with the argument in RefJ^ 
The direct thermal force would become physical, if we 
consider an ensemble of domain walls which are allowed 
to annihilate and to be created (see Sec. IVIIII) . 

The thermal effect for a domain wall therefore arises 
only from the magnon interaction, (Eq. dT51) l. Not- 
ing that ^ sin0Vi0, the force and torque induced 

by the magnon interaction read 


F„ 


SHf 

SX 

SHf 

5(t) 


= 0 

2JS . 

Jm.. 



A 


2N^JSa^ . 
A 

(29) 


In Eq. (l29l) . the force from the magnon disappears since 
we consider here a static one-dimensional domain wall, 
in which case the magnons are in the perfect adiabatic 
limit^; In other words, magnon modes are orthogonal 
to the collective modes of a static one-dimensional do¬ 
main wallas. The force from the magnons may arise when 
the wall is dynamic^ or higher-dimensional and in the 
presence of dipolar interaction^. Here we introduce phe¬ 
nomenologically a thermal force by use of a dimension¬ 
less parameter /3t as Fp = zT. Including 

this force and the thermal magnon effects, the equation 
of motion for the wall is 


where tan ^ and c is an integer representing the chi¬ 
rality. Using the collective coordinates, X{t) and U(t), 
representing the center of the vortex in two-dimensions, 
the spin structure is S{r,t) = Sv{r — X{t)), where 
X{t) = {X{t),Y{t)). Similarly to the domain wall case, 
we consider the case where vortex cannot annihilate; All 
the thermal effect then arises form magnons. The spin- 
transfer torque due to magnons, Eq. m , is written 
when magnon current is uniform as 

where 6S is the change of S when the origin of X is 
shifted by an amount SX = {SX,6Y). (This represen¬ 
tation is convenient to focus on physical contribution in 
topological quantities for spin^. ) The magnon current 
density is the two-dimensional one = jmLz 

where is the thickness of the system). The force due 
to the spin-transfer is therefore 


_ SHf 2J 

SX, S'2 


d\S ■ [X,S X 


(34) 


This force is characterized by the topological number 
density of the vortex, 


h f d?r 


2TrhS 


G = ^ I ^S-[X,.SxXyS]= , 


(35) 


as 


TT. — _ ri X 

^ ^ ^ Jm ’ 


(36) 


i + a- 
X — aX(j) 




hX 


K±XS 

2h 

K±XS 

2h 


sin 2(j) + 
sin 2(j) — 


X 

hN^S^'^ 

2Ja^ . 


(30) 


The magnitude of the magnon current is evaluated by 
use of a linear response theory with respect to FIat 
Sec. m 


VIII. VORTEX 


Let us next consider a case of a single vortex which 
appears in a two-dimensional sub micron size disks. A 
Hamiltonian describing the vortex is 


H = 





(31) 


and the vortex structure with vortex number of unity is 
approximately represented as^ 


where G = Gz. The equation of motion of a vortex is 
then^i^ 


-G X X + aDX ^ (37) 


where D cx ^ j cPr sit? 9{V fpY' is a form factor for the 
damping. The equation is written as 


xYJf 

ri 

The velocity is thus 
1 


aD 


= -^G X X. 


X = - 




aD\ 2Ja? 


G 


(2d) 


(38) 


(39) 


where ay is the y-component of the Pauli matrix. We 
consider the case of temperature gradient along the x- 
direction. Since the magnon current is along VT, the 
vortex velocity is 


A 

Y 


1 2Ja? ,( 2 d) 



(40) 


Sy{r) = S ^cos e^; -|- sin , (32) 


In the study by a numerical simulation and the Eocker- 
Planck equation in Refj^i carried out for a skyrmion, 
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there was a direct force term proportional to Va;T. (A 
single skyrmion and a vortex are described the the same 
equation called the Thiele equation^ and thus have the 
same response to applied forces.) Such direct force arises 
in our formalism by taking account of annihilation of vor¬ 
tices and by regarding the energy current, Eq. da , as 
physical one contributing to the entropy production. The 
energy current density in this case becomes js = 7 vA' 
( 7 v is a temperature-dependent coefficient), and the di¬ 
rect thermal force becomes 

Ft = = ^VT. (41) 


The lessor Green function is written in terms of retarded 
and advanced Green’s functions as^ 


< 


9q,u:-^9q^u:+^ 


— n 



g 

^ q^uj — 


^9 


2 


a 


2 


— n 



9 


r 

q,u) — 


n9 


2 


r 


n 

2 



9 


a 

q,uj — 


0.9 


2 


a 

q 


,cj+ 


n 

2 




^ '( 

i^) [9 


a 

q,iij 



Since the annihilation of vortices requires a finite excita¬ 
tion energy, vanishes quickly at T = 0 , and the force 
would vanish at T = 0. It would be interesting to study 
whether the direct thermal force exists or not experimen¬ 
tally. 


where n{uj) = — 1]“^ is the Bose distribution func¬ 

tion, n'{u}) = ^, /3 = l/{kBT) {ks is the Boltzmann 
constant) and we have neglected contribution of the or¬ 
der of The retarded and advanced Green’s function 
are 


IX. LINEAR RESPONSE THEORY OF 
MAGNON CURRENT 


9quj 


1 

LO — UIq + iauj ’ 


(46) 


We calculate the magnon current induced by the tem¬ 
perature gradient within the linear response theory. We 
carry out the calculation by use of the non-equilibrium 
Green’s function. (Standard Kubo formula calcula¬ 
tion is also applicable.) The interaction Hat in the 
Fourier representation of the magnon operators {b{r, t) = 



(42) 


and gq^ = (gq^)*, where uiq is the angular frequency of 
magnon with a wave vector q and the effect of Gilbert 
damping is included as an imaginary part. The ‘param¬ 
agnetic’ magnon current is thus 

« + * ^ AT,*(-r!)wn(cu) , 

2 ’ q 

(47) 


where fl is an infinitesimal external angular frequency 
and V is the system volume. The ‘paramagnetic’ and 
‘diamagnetic’ magnon currents, Eq. m and (E21), re¬ 
spectively, are 



(43) 


where we used 2JSqj{gq^^)‘^ = dq^gq^. The ‘diamagnetic’ 
magnon current is similarly calculated as 



(48) 


We see that the ‘diamagnetic’ current cancels an equi¬ 
librium contribution of the ‘paramagnetic’ current, re¬ 
sulting in (assuming rotational symmetry and using 
f ^(-ifl)AT,i(-^2) = At = ^) 


where = -i (h\,^,hq^ 


is the lesser component 

•fi in 1■I'-vn rvnCi53 


ear order in Hah 


as 


.( 0 ) 

1 = 


2JS 


f duj f dft ^ 
q 


9q,u>-%9, 


^q,uj+^ 


Jm,z — iT, 

where the coefficient k is 




lin- 

JS 1 

K =- 

m T 

f duj 1 

j 

- < 

AJS 

2 1 1 

(44) 

3h 

“ TV 


(49) 


[{u: - ujqf + {aujff' 
(50) 
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To evaulate the summation over q, we consider the 
case of uniaxial anisotropy {K_i_ <C K) for simplicity . 
The magnon energy is then hiOq = JSq^ + Agw, where 
Asw = KS is the gap of spin wave. The summation 
over q is carried out in three dimensions. Defining e = 
JSq^ + Asw, we obtain 


1 




1 


212 




1 


47r2(J5)5/2 

1 

2-k{J 


poo 

/ de(e-Asw)^^^ 

Asw 

e(u: - Asw) 


[(e - w )2 + (aw)2]2 
(w - Asw) 3/2 


(51) 


where 0(x) is step function and we used an approximation 
(using a <C 1) 

1 27r 

77-7TT7—7Ti2 - 7—(^2) 

[(e — w)2 + (Q;a;)2J2 [aujy 

The result is thus 


1 1 
'37r2(J5)3/2^ 


dniio) 

duj 


27r2(J5')3/2 a 


fcB(fcBT)i/2i^(/3Asw), 


(53) 


where 


poo 

F(^Asw)= / 

•1/3 As 


d Jx /lAsw)^/^ ^ 

- 1 


At high temperatures, /lAgw >C 1, A = 2.315. We 
see that the thermal transport coefficients are propor¬ 
tional to the damping time constant, o;”^, in the same 
manner as in the electric transport coefficients in met¬ 
als are proportional to the elastic lifetime, r. The result, 
K oc a/T/ a, agrees with a semi-classical calculation in the 
supplement of Ref JA. 


X. DOMAIN WALL SOLUTION 


Using the result of magnon current, Eq. with Eq. 
(15311 . the equation of motion for a domain wall, Eq. dsni), 
reads 


where 


and 


1 , ^ Q '^T 

p + a— = Pt — 

(55) 

X — aXp = Vc sin 2p — FtUt, 

(56) 

KtXS 

(57) 

“ 2/1 ’ 

n 

(58) 


represents the scale of the velocity induced by the tem¬ 
perature gradient (with positive direction chosen as the 
direction from the high- to the low-temperature region). 
The spin-polarization coefficient in the thermal magnon 
spin-transfer is temperature-dependent as 


_ 2Jk 
kso^ 


F 1 

7T^\/S Oi 



(59) 


where F = 2.315 considering high temperature (/3Asw “C 
1). Note that there is a minus sign in the magnon spin- 
transfer term (proportional to Ft) in Eq. (l30l) indi¬ 
cating that the magnon spin-transfer pushes the wall 
to the hotter end of the system as has been noted 
previouslyi^ii^. Equation (15^1) indicates a clear distinc¬ 
tion between the magnon spin-transfer effect from the 
electron spin-transfer effect in metals, namely, the spin- 
transfer efficiency parameter. Ft, grows at small damp¬ 
ing (a), while the parameter is independent on a in the 
case of electron spin-transfer effect. 

The equations of motion, Eq. dsni), has the same form 
as the current-driven case in metals. As was shown in 
Refii^, the behavior of the wall are distinct at small and 
large ut compared to the crossover velocity (the Walker’s 
breakdown velocity), defined as 


Ft + ^'' 

' rv 


(60) 


which depends on the temperature and a. At small driv¬ 
ing velocity, IutI < |uc|, the tilting angle of the wall, 0, 
reaches the terminal angle determined by 

sin2(() = ^ , (61) 

V a J Vc 


and the terminal velocity of the wall is a constant. 


Uw 



a 


(62) 


The wall thus moves from the high- to low-temperature 
region in this case, due to a pressure on the wall. Above 
the crossover velocity, \ut\ > |uc|, the angle cj) becomes 
time-dependent and the average wall velocity is 


Uw 



1. (63) 


Eor UT » |uc|, Ww = i^{-FT+al3T)uT- If Ft- ajdr > 
0, the direction of domain wall motion therefore changes 
around ut ~ 0(|mc|); at small VT, the motion is towards 
the colder end while it is opposite in large VT region, 
dominated by magnon spin-transfer effect. At ut ^ |uc|, 
the wall velocity is Uw/VT = —^^^{Ft — a/Sr)- 

The domain wall velocity as a function of ut/vc is plot¬ 
ted in Fig. [1] for different temperatures and Pt- When 
there is no force {Pt = 0), there is a threshold value of 
for Ut for thermal motion to set in, and the motion 
























FIG. 1. (a-c): The domain wall velocity in the unit of 

Vc as a function of ut/vc for different temperatures and for 
pT = 0, 0.005 and 0.02 at a = 0.0075. Negative v„ is towards 
the direction to the hotter end. Parameters used are a = 12A, 
S = 14, = 4.2 X 10^ K. The case of a = 0.02 with 

Pt = 0.005 is plotted as (d). Effects of extrinsic pinning 
is not taken into account. It is seen that the magnon spin- 
transfer effect is enhanced at high temperature and when a 
is small. 



FIG. 2. The domain wall velocity for Pt = 0 and P = 0.005 
at T = lOOK and 300K at a = 0.0075 for (a) small and (b) 
large utIv^. We see that dependence on value of Pt is weak. 


towards the hotter regime occurs above the threshold. 
The wall speed increases as function of the temperature, 
since the effective spin polarization, Pt, increases at high 
temperatures rather significantly (oc VT). If Pt is pos¬ 
itive, wall motion towards the colder regime occurs at 
small ut{< ImcDj and the magnon spin-transfer regime 
appears at IutI |ucl- For large a, the wall speed is 
smaller (Fig. (Hd)). 

In YIG, a = 0.0075, a = 12A, and the magnetization is 
M = 1.4 X 10® A/m. If we convert the magnetization to 
the average magnitude of spin by M = we obtain 

the effective spin of ^ = 14. The exchange stiffness is A = 
4.7 X 10“^^ J/m, which corresponds to = Aa = 5.8 x 
10“^^ J = 4.2 X 10^ K). In the three-dimensional 

case, we therefore have Pt = 0.40 x y^r(K). (At T = 300 
K, Pt = 6.7 and the crossover velocity is Uc = —Vc x 
0.15.) The absolute value of the wall speed depends on 
the hard-axis anisotropy energy, K± (vc oc K±). The 
easy-axis anisotropy energy of YIG is KS^aP = 1.2 x 
10“^^ J, and the hard-axis energy is much weaker than 
the easy-axis one. Let us choose K±/K = 10“®, i.e., 
Vc = 2.4 x 10“^ m/s. We consider the case /3 t = 0 first. 
We have |mc| = 0.36 x 10“^ m/s at T = 300 K. In this 
case, a temperature gradient of VT = 20 K/mm in Ref.— 
corresponding to ut = 3.8 x 10“® m/s (ut/w = 0.16) is 
close to the threshold {ut — |uc|), and the wall speed is 
read from Fig. [T](a) and Fig. ^a) as Ww = —0.37 tc = 


—0.89 X 10“^ m/s. In the case of Pt = 0.005, ut/uc = 
0.16 leads to = —0.50z;c = —1.2 x 10“^ m/s (Fig. [TKb) 
and Fig. Ha) ), and we see that the order of magnitude 
of the speed does not depend strongly on Pt (Fig. Ha))- 
An interesting fact is that the wall velocity changes sign 
as we lower the temperature if Pt is finite. In fact, when 
we set utJvc = 0.16 for Pt = 0.005, the wall velocity is 
negative at T = 300 K, while it is positive at T = 100 
K as seen in Fig. Ha)- The sign change occurs around 
T = 190 K, where Uc{T) ~ ut- 

If hard-axis anisotropy is weaker, K^/K = 10“^, Uw = 

— IlUc = —2.6x 10“^ m/s at r = 300 K and = —6uc = 

— 1.4 X 10“^ m/s at T = 100 K for ut = 3.8 x 10“® m/s 
and for both /3t = 0 and 0.005 (Fig. H^))- 

In the experiment on YIG at the room temperature in 
Ref—, the wall moved to the hotter regime at speed of 
1.8 X 10“^ m/s at VT = 22 K/mm. This speed is smaller 
than our theoretical estimate by two orders of magnitude, 
and we suspect that the small value would be due to the 
extrinsic pinning. In fact, the experiment was carried 
out under rather strong field of R = 6 x 10“® T (60 Oe), 
which adds a large additional force term = 1 x 10® 
s“^ to the right-hand side of Eq. (ESI). This force is 
expected to have been necessary to compensate extrinsic 
pinning effects, and this fact and the avalanche behavior 
of walls observed suggest a strong influence of extrinsic 
pinning in their sample. 

In the experiment carried out in metals in Ref.—, the 
applied VT was order of 100 K/mm. The motion was ob¬ 
served in the presence of ac magnetic field to remove the 
extrinsic pinning effects and wall velocity of 2 mm/s to¬ 
wards the colder end was obtained. The direction would 
be understood by the spin-transfer effect of thermally- 
induced conduction electron flow. 


XI. CONCLUSION 

To summarize, we applied a thermal vector potential 
theory developed in Ref— to describe magnetization dy¬ 
namics in ferromagnetic insulator driven by a temper¬ 
ature gradient. We have evaluated the magnon current 
induced by the temperature gradient within the linear re¬ 
sponse theory and it was found to be proportional to the 
inverse of the Gilbert damping parameter. The effect of 
magnon current is thus dominant in weak damping sys¬ 
tems as was argued in Ref—. The magnon current was 
shown to exerts a spin-transfer torque on a domain wall, 
and this effect tends to drive the structure towards the 
hotter side. The case of a vortex (or a skyrmion) was 
also discussed. 
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Appendix A: Energy current density of localized 
spin 


where we used the fact that spins on different sites com¬ 
mute each other. The contribution from the last two 
terms vanish after summation over a and cr', and we ob¬ 
tain 

cta'fS/S'S cr,cr'=± 

^ i^i+<7l3+a'l3' ~ ^i+a'P')- (^6) 


We derive here the expression for the energy current 
density of the localized spin by evaluating the time- 
derivative of the energy density quantum mechanically. 
The Hamiltonian we consider is the one with the ex¬ 
change interaction and easy and hard axis anisotropy 
energies, Eq. We carry out the calculation on a dis¬ 
cretized lattice, since the estimation of commutators goes 
straightforwardly, and then take the continuum limit. 
The energy densities at site i are 


= 

S^ = 


^ \ ' \ ' QOL na 

2^5 ^i+crP 


OL0 cr=± 




£ 


K, 




(Al) 


where a, /3 runs x, y and z, and a denotes the positive and 
negative directions. Corresponding terms of discretized 
Hamiltonian are 


The commutators including the exchange interaction and 
easy-axis anisotropy energy are 




4a5 


OCX OCX QZ 

*^z+(T^‘^2+cr/3 


cxcx' j3 <7 — -^ 


(A7) 


The hard-axis anisotropy contribution has the same form 
but with the index 2 ; replaced by y. 

Let us take the continuum limit. Contribution from 
the exchange interaction is {y. = x, y, z) 

act.'^5 fT=± 

= -tnv-ji\ (A8) 




K 
i 5 


i 


t 


* (A2) 


We hrst study the exchange interaction contribution 
by evaluating 

j aa'j3^' aa' 

(A3) 

By use of [S°‘,S^] = iSij€afi^Sj and 


where 


jii = 


ha^ 


[(V^S) X S] 


(A9) 


is the energy current from the exchange interaction. Sim¬ 
ilarly 


[£;',HK] + [£!^,Hj] = -iHy-ji^ 

[£/ ,HkA + [ff", Hj] = -thV- jE- , (AlO) 


where 


[AB, CD] = A[B, C]D + [A, C]BD + CA[B, D] + C[A, D]B, 

(A4) 

we see that 

E r QQ! Qa Qa' oa' ] 

j 

6 

I QC>t' QC>t qS I Qa' QC>t \ 

— ‘^2+cr/3 — ^i+alS j 

= I E 

5 ^ 


(A5) 


■JK ^ JK 
— 


{SxV^S-y^Sx S),S, 


Aha^ 

S^{S X V^S-V^S X S), 


— 


JK 


[S X X S)ySy 


4fta3 

+ Sy{S X S)y 


(All) 


In the classical case, where ordering does not matter, we 
obtain 


JK, 


J . 


38^^ = X = -7g^V^5 • {Bk X S) 


J , 


3e!;^=-J,^^,S-{BK^xS), 


(A12) 
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where 7 g is gyromagnetic ratio and 

= -KS^z 

= K^SyV, (A13) 

are the effective magnetic fields for the easy and hard axis 
anisotropy energies, respectively ( y and z are the unit 
vectors along y and 2 axis, respectively). There is no con¬ 
tribution to the energy current from the commutators of 
anisotropy energies, since local quantities commute each 
other. The total energy current for the Hamiltonian ([S]) 
is therefore 

js,ii = + 3s^ + 

= -Je^y^S-{BHX S), (A14) 


where 

S TT 

n-t^BH = j^=-JV^S-KS,z + K^Syy, (A15) 

is the total effective magnetic field for the Hamiltonian 
H. The equation of motion for S is 

S = 'ygBx X S, (A16) 

and thus we finally obtain the energy current density as 

js,^. = -^V^S-S. (A17) 

This form in fact the one obtained from a general def¬ 
inition in terms of the Lagrangian L, Tqi = (dtS)-^^. 
It also agrees with the one proposed in Refi^ based on 
a symmetry argument, although the spin relaxation was 
believed there to be essential for emergence of a term of 
Eq. (IATtIi . 
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